Thermodynamic properties and shear viscosity over entropy density ratio of nuclear 
fireball in a quantum-molecular dynamics model 
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Thermodynamic and transport properties of nuclear fireball created in the central region of heavy- 
ion collisions below 200 MeV/nucleon are investigated within the isospin-dependent quantum molec- 
ular dynamic (IQMD) model. These properties include time evolutions of the density, temperature, 
chemical potential, entropy density (s) and shear viscosity (rf) as well as density and temperature 
dependencies of the ratio of shear viscosity over entropy density (j]/s) etc. Based on the shear 
viscosity parametrization developed by Danilewicz and entropy density which is obtained by a gen- 
eralized hot Thomas Fermi formalism, the ratio of shear viscosity over entropy density is calculated 
in the whole collision process as well as in the freeze-out stage. With the collision goes on, a tran- 
sient minimal rj/s with the value around 5/4-k occurs in the largest compression stage. While, the 
relationship of rj/s to tempertaure (T) in the freeze-out stage displays a local minimum which is 
about 9-10 times l/4n around T = 9-15 MeV, which can be argued as indicative of a liquid gas 
phase transition. In addition, the influences of nucleon-nucleon cross section (onn) and symmetry 
energy coefficient (C s ) are also discussed, and it is found that the results are sensitive to cmn but 
not to Cs- 
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I. INTRODUCTION 

In the past decades, extensive experimental and the- 
oretical efforts have been devoted to search for the nu- 
clear liquid-gas phase transition (LGPT) in intermediate 
energy heavy-ion collisions (HIC) jl|-fll|. Many probes 
have been suggested for the onset of nuclear LGPT. For 
instance, the fragment size distribution [12| and its rank 
distribution [ll[, the largest fluctuation of the heaviest 
fragment fl3j , caloric curve dH, bimodality [l| etc. In 
addition, it has been observed that the ratio of shear vis- 
cosity to entropy density (rj/s) reaches its local minimum 
at the critical temperature for a wide class of systems. 
For instance, empirical observation of the temperature 
or incident energy dependence of the shear viscosity to 
entropy density ratio for H2O, He and Ne2 exhibits a 
minimum in the vicinity of the critical point for phase 
transition [l5[. And a lower bound of 77/s > 1/Att ob- 
tained by Kovtun-Son-Starinets (KSS) for infinitely cou- 
pled super-symmetric Yang-Mills gauge theory based on 
the AdS/CFT duality conjecture, is speculated to be 
valid universally 1|| In ultra- relativistic HIC 0- 
[23j . people have used the ratio of shear viscosity to en- 



tropy density to study the quark-gluon plasma phase and 
the extracted value of 77/s seems close to the KSS bound 

(1/47T). 

So far there are many interesting investigations on the 
ratio of 77/s, but it is still rare to study the behavior 
of 77/s during the heavy- ion collision at intermediate en- 
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ergies. Furthermore, the influences of nucleon-nucleon 
cross section and nuclear symmetry energy are less dis- 
cussed. 



In this work we use a microscopic transport model 
known as the isospin-dependent quantum dynamics 
model [H[ to simulate Au+Au central collisions. In order 
to study the effect of nucleon-nucleon cross section, 0.5 
times and 1.5 times normal nucleon-nucleon cross section 
are also used in the simulation. On the other hand for 
the symmetry energy which is importan for asymmet- 
ric nuclear matter as well as nuclear astrophysics [34l - 
[36j . different symmetry energy parameter is employed, 
namely 15 MeV, 25 MeV and 35 MeV. The generalized 
hot Thomas Fermi formalism (GHTFF) ^MX^ and the 
transport formula [25| are employed , respectively, to ex- 
tract thermodynamic and transport properties of the nu- 
clear fireball which is located in the central region with 
a moderate volume. At the end, different correlations 
between the extracted thermal and transport properties 
are discussed and we find a good agreement with our 
previous calculations (30j . 



The paper is organized as follows. Section 2 provides 
brief introductions for the IQMD model, GHTFF as well 
as transport formula for shear viscosity. In Section 3 we 
present the calculation results and discussions, where the 
time evolution of thermodynamic quantities and shear 
viscosity over entropy density are focused. Fianlly a sum- 
mary and outlook is given. 
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II. 



MODEL AND FORMULISM 



A. Quantum Molecular Dynamics Model 

The quantum molecular dynamics (QMD) [52[ 
model approach is a many-body theory which describes 
heavy ion collisions from intermediate to relativistic en- 
ergy. The isospin-dependent quantum molecular (IQMD) 
[4j],|42[ model is based on the QMD model, including the 
isospin effects and Pauli blocking. Each nucleon in the 
colliding system is described as a Gaussian wave packet, 



jexpl-p^t) -r. 



(2ttL) 3 / 4 ~"~" l 7 
4Z 



(1) 



Here Ti(t) and Pi(t) are the mean position and mean 
momentum, and the Gaussian width has the fixed value 
L = 2.16/m 2 for Au + Au system. The centers of these 
Gaussian wave packets propagate in coordinate (R) and 
momentum (P) space according to the classical equations 
of motion: 



0(H) 



0(H) 
dpi 



(2) 



where (H) is the Hamiltonian of the system. 

The Wigner distribution function for a single nucleon 
density in phase space is given by 



/i(r,p,0 



exp 



■ exp 
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(3) 



The mean field in IQMD model is written as 

U (p) = U Sky + J/coui + U Yuk + U sym , (4) 

where f/skyi Uqou\, Uy u k, and U sym represents the Skyrme 
potential, the Coulomb potential, the Yukawa potential 
and the symmetry potential interaction, respectively (4[J . 
The Skyrme potential is 



U S ky = a(p/po) + /3(p/po) 7 , 



(5) 



where po = 0-16 fm -3 and p is the nuclear density. In the 
present work, the parameters a = —356 MeV, j3 = 303 
MeV, and 7 = 7/6, correspond to a soft EOS, are used. 
U Yuk is a long-range interaction (surface) potential, and 
takes the following form 

jjYuk = (y y /2)^exp{Lm 2 )/r rj 

\exp{mrij)erfc{^/Lm — rij/\AL) 
— exp(mrij)erfc(y/Lrn + rij/V^L)], (6) 



1.25fm-\ L = 2.16 fm 2 , and 



with V y = 0.0074GeV, m 

is the relative distance between two nucleons. The 



symmetry potential is U s , 



C K 



-Pp 



where p n , p p , 



' s V m ~» po 

and po are the neutron, proton and nucleon densities, 
respectively. C s is the symmetry energy coefficient, and 
three different values of 15, 25 and 35 MeV are used to 
study the influence on the ratio of r\js. 

Furthermore the isospin degree has entered into the 
cross sections, and is similar to the parametrizations of 
VerWest and Arndt, see Ref. (HT| ■ The cross section for 
the neutron- neutron collisions are assumed to be equal to 
the proton-proton cross sections. In order to study the 
effect of cross section on the ratio of rj/s, the nucleon- 
nucleon cross section is multiplied by a parameter C a . 
Three different situation is considered when C a equals 
0.5, 1.0 and 1.5 respectively . 



free 
NN ■ 



(7) 



From Eq. (3) one obtains the matter density of coor- 
dinate space by the sum over all the nucleons, namely 

At+Ap 
0=1 

A T +Ap 



1 (r-rj(f)) 2 



2L 



The kinetic energy density in coordinates space could 
also be calculated from Eq.(4) by 



PK(r,t) = ^ 



2m 



Pj{r,t). 



(9) 



B. The Generalized Hot Thomas- Fermi Formalism 

The thcrmodynamical properties of hot nuclear mat- 
ter formed in heavy ion collisions, e.g. temperature and 
entropy density, can be extracted by using the approach 
developed by Faessler and collaborators [37H391 l43l |44| . 
In this approach one starts from a microscopic picture 
of two interpenetrating pieces of nuclear matter and de- 
duces the thermal quantities from the matter density and 
kinetic energy density obtained during the collisions. In 
this paper, the extraction of thermal properties of the 
hot nuclear matter is done in two steps. First, based 
on the IQMD simulation, one could calculate the nuclear 
matter and kinetic energy densities at each point in co- 
ordinate space at every time step. Second, by employing 
the hot Thomas-Fermi formalism, we could obtain the 
corresponding thermal properties for every set of nuclear 
matter density and nuclear kinetic energy density. De- 
tailed information can be found in Refs. |37l. l38j. 

In GHTFF, the local nuclear matter density pi is ex- 
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pressed by 



1 . , 1 ,2mT x o /? 
x[f( (H ,K 0i ) + Ji /2 (jM,K^)], (10) 



where p,i is the local chemical potential and T is nuclear 
temperature. The function / and K will be given later. 
The local kinetic energy density reads 



2m 

1 ,2mT, 



TifjH) = ^ {p l ) + —{ ^ 2 



5/2 



x [l^oifilH, K m ) + -J 1/2 {fXi,K m ) 



K 0i 



Ji(fj,i - x 2 )dx] + An(/ii). (11) 
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And the entropy density s» is written as 



/ \ 1 / \ 1 ,27TLT \q/9 

Si (Mi) = 2 S o(^) + ^2(~p-) 



-MiJ 1/2 (Mi,-^oi) 



2 / Ji(^ ? ; -cc 2 )d.T]. (12) 
/o 



Here i = 1,2 represents the projectile and target, and 



A". 



01 



i , A 02 = K R - A i with K r 



TiKr. 



/2mT' " u± ' \/2mT' 

Jn{p>) = Jn(fJ-, oo) is the Fermi integrals, i.e. 



and 



iq 1 + cxp(x — /i)' 
and /(/ii, A 0i ) = A ^n[l + exp(/^ - Ati(#i) = 0, 



FIG. 1: (Color online) Shear viscosity of nuclear matter as a func- 
tion of p/po and T with the Eq. l|13|l . Different colors represent 
different p/po, which arc illustrated in the inset. 



C. Shear Viscosity Formalism 

For largely equilibrated systems, fluxes of macro quan- 
tities, leading to dissipation, are proportional to gradi- 
ents within the system. The shear viscosity denoted as 
r\ is the coefficient of proportionality between anisotropy 
of momentum-flux tensor, including dissipation and ve- 
locity gradients. In the Boltzmann statistical limit the 
shear viscosity corresponds to the first order Chapman- 
Enskog coefficients. In Ref. [25|,[2(| the nuclear shear vis- 
cosity has been derived from the microscopic Boltzmann- 
Uehling-Ulcnbcck equation and can be parameterized as 
a function of density p and temperature T: 



At 2 (/U 2 



x[Ji(/i 2 ) - JiGu 2 , A 2 2 ) - A 02 /(ii 2 , A 02 )] 
+k%p 2 (fj, 2 ), 



where k R represents the relative momentum between the 
target component and projectile component. From the 
cq. (|10lllll2|) . one can obtain the thermal properties by 
inversion in principle. But such an inversion procedure 
is practically not feasible due to the complexity of the 
equations. Therefore, a more practical way is chosen to 
solve this problem. First we calculate all reasonable com- 
binations of temperature (T), relative momentum (Kr) 
and chemical potential (pi), and make a table of all the 
related values (T,Kn,/j,i,pi,Ti,Si). Second from the nu- 
clear matter density and kinetic energy density eq. (|8!9l) . 
which is extracted from the IQMD simulation, we can get 
the thermal properties by comparing the calculated ta- 
ble. More detailed information can be found in Ref. 1371. 



V(~,T) 
Po 



T 2 pa 

5.&VT 

+ 1 + 160T 



22 / P nO.7 

T 2 10- 3 >o 



-2 ' 



(13) 



-3 



where r\ is in MeV/fm 2 c, T in MeV, and po=0-168fm 

Figure Q] shows rj as a function of T and p/po- One can 
see that rj exhibits a very distinct minimum when nu- 
clear matter density is less than normal nuclear density. 
And as the density increases, the critical temperature 
also get larger, e.g. for normal density the critical tem- 
perature locates around 10 MeV, but for 4 times normal 
density it is almost 80 MeV. This conclusion is coincident 
with macroscopic result. For the equilibrium of the con- 
sidered nuclear matter at the very starting stage being 
poor, the shear viscosity extracted should be treated as 
the transport properties of largely equilibrated nuclear 
matter with the same density and kinetic energy density. 
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III. CALCULATION AND DISCUSSION 

In present work, we simulate head-on collision of Au 
+ Au at different beam energies. The reason why we 
choose the central collision is that the participant zone is 
the maximal and the matter is more nuclear liquid-like 
during the early time evolution of collision. 

A. Time evolution of density profile in reaction 

plane 

Fig. [2] shows the evolution of density profile in X-Z 
plane for Au+Au head-on collisions at 50 MeV/nucleon 
and 200 MeV/nucleon at 2, 10, 20, 30, 50 and 80 fm/c. 
The zero point of time is set at the initial contact between 
project and target (the first left panel) . With the collision 
goes on, the system overlaps and seems more isotropic in 
phase space. In order to calculate the thermodynamic 
quantities in different time steps, we select the central 
sphere with radius r=5 fm, which is the volume of nuclear 
fireball in this paper. 

B. Time evolution of thermodynamic variables 

The time evolution of the average nuclear matter den- 
sity (panel (a)) and kinetic energy density (panel (b)) in 
the central region, with C a = 1.0, is showed in Fig. [31 It 
is interesting that the maximum density reached is about 
1.5/30 to 2.0po an d the maximum kinetic energy density is 
10 McV fm -3 to 25 McV fm -3 for the energy displayed in 
the picture. Along the time scale of the collision one can 
see that both p/ po and r are reaching their maxima at 
about 20 fm/c and a bit earlier for higher energy. After 
the compression stage the matter starts to expand and 
flow out of the central region, mainly in the transverse 
plane, the matter density drops to very small values and 
the central region is cooled down. But the hot and dense 
nuclear matter exists much longer when the incident en- 
ergy is low. At about 80 fm/c the hot and dense matter 
disappears. 

Time evolution of matter density and kinetic energy 
density are shown in Fig. @J when C a and C s are set 
by different values . The upper panels demonstrate the 
different cross section situation (C a € [0.5,1.0,1.5]), on 
the other hand the bottom panels are for different sym- 
metry energy. In panels (a) and (b), it is easy to find 
that the nuclear matter density and kinetic energy den- 
sity are different to each other when C a is different. We 
found that there is no difference for the nuclear matter 
density during compression stage. But when the system 
starts to expand, the larger nucleon-nucleon cross section 
makes the dense matter stay longer. In contrast with the 
behavior of density around the maximum compression 
stage, more distinction for the kinetic energy density is 
exhibited. It shows that the smaller the nucleon-nucleon 
cross section, the larger the kinetic energy density. But 



these curves almost overlap each other after 40 fm/c, it 
may be understood that the longitudinal energy enters 
the central region more easily when the nucleon-nucleon 
cross section is small. It should be noted that the large 
kinetic energy density dose not mean higher tempera- 
ture, since the kinetic energy is not calculated in the 
center of mass frame, detailed information can be found 
in Rcfs. [13-111, SSI- In addition, the extracted den- 
sity and kinetic energy density show insensitivity to the 
symmetry energy as depicted in panel (c) and (d). The 
curves are overlapped to each other for the whole process. 
As has been discussed in the previous paragraph, this 
leads the nuclear matter to the same thermal properties. 
So the extracted thermal properties based upon the hot 
Thomas-Fermi formulism is exactly the same with each 
other. So in this paper we just investigate the nucleon- 
nucleon cross section effect on thermodynamic and trans- 
port quantities in the following texts. 

The time evolution of temperature is plotted in Fig. [5] 
Fig. [SJa) shows the time evolution of temperature with 
(7,7=1.0 at different incident energies. Fig.^b) shows the 
time evolution of temperature when the incident energy 
is fixed at 130 MeV/nucleon but with different cross sec- 
tions. The following pictures are arranged with the same 
mode, i.e panel (a) represents a constant cross section 
C CT= i.o but at different incident energies; panel (b) means 
a constant incident energy at 130 MeV/nucleon with dif- 
ferent cross section. For a given beam energy, temper- 
ature increases at first, then reaches a local maximum 
about 20 fm/c and decreases till a saturate value about 
80 fm/c. The higher the incident energy, the larger the 
maximum value. The corresponding time at maximum 
value is a little earlier than that for the density and ki- 
netic density. In panel (b), it is found that the larger 
cross section makes the system a little hotter. The rea- 
son is that there are more frequent nucleon-nucleon col- 
lisions as (Jnn becomes larger, which makes the transla- 
tion from the longitudinal energy to thermal energy more 
efficiently. 

Fig. [5] shows the time evolution of chemical potential 
(p). Again, the left panel displays the normal cross sec- 
tion one, we can find that p increases in the compres- 
sion stage and decreases in the expansion stage, and the 
lower the incident energy, the larger the chemical po- 
tential. This might be understood as a large compound 
nucleus is formed during the compression stage, and the 
lower the incident energy, the larger the compound nu- 
clei is. In panel (b), it shows that the chemical potential 
becomes generally larger when nucleon-nucleon cross sec- 
tion is larger. 

The time evolution of entropy density is plotted in 
Fig. [7] It is found that the entropy density almost syn- 
chronically evolves with the temperature. The higher 
the incident energy and nucleon-nucleon cross section, 
the larger the entropy density is. 
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FIG. 2: (Color online) X-Z density profile in different time step for Au+Au head-on collisions at 50 McV/nucleon (upper panels) and 200 
MeV/nucleon (lower panels), respectively. From left to right panels, time step is 2, 10, 20, 30, 50 and 80 fm/c, respectively. 



C. Ratio of shear viscosity to entropy density 



Now we can move to the discussion on transport co- 
efficient. Since the nuclear participant in central region 
could be seen as nuclear fluid, we adopt the Eq. [T3] to 
calculate the shear viscosity. Unlike the Green-Kubo 
formula [3l|, HH, the advantage of the equation is that 
we can investigate the time evolution of shear viscosity 
in the framework of transport model. But it should be 
noted that Eq. Q2]is principally applicable when the sys- 
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FIG. 3: (Color online) Time evolution of mean matter density (a) 
and kinetic energy density (b) at different beam energies. 
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FIG. 4: (Color online) Time evolution of mean nuclear matter 
density ( (a) and (c)), and the kinetic energy density ((b) and (d)) 
in the central region defined as a sphere with radius equals 5 fm of 
head-on Au+Au collisions at 130 MeV/nucleon. Different nuclcon- 
nucleon cross section ((a) and (b)) and symmetry energy parameter 
((c) and (d)) are used. 



tern is largely equilibrated. However, a full equilibrium 
is hardly achieved during the whole heavy-ion process. 
So the shear viscosity extracted here should be treated 
as the properties of an equilibrated nuclear fireball with 
the same thermodynamic state as the simulated one. 

Fig.[5]displays the time evolution of shear viscosity (77), 
it shows an increase in earlier stage and then drops with 
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FIG. 5: (Color online) Time evolution of temperature generated in 
the central region at normal nucleon-nucleon cross section at differ- 
ent inciden energies (a), or at 130 MeV/nucleon but with different 
&NN (b). The incident energies and NN cross section parameters 
arc illustrated in the inset. 
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FIG. 6: (Color online) Same as Fig. [5] but for chemical potential. 
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FIG. 7: (Color online) Same as Fig. [5] but for entropy density. 



FIG. 8: (Color online) Same as Fig. [5] but for shear viscosity r\ of 
the central fireball. 



time. As Eq. [JJ shows, here the shear viscosity depends 
on both temperature and density which varies with time. 
Roughly speaking, the shear viscosity increases in the 
compression stage and decreases as the system expands. 
The smaller the nucleon-nucleon cross section, the larger 
the viscosity in the maximum compression stage. 

When the entropy density is taken into account, the ra- 
tio of shear viscosity to entropy density shows a minimum 
near maximum compression point as shown in Fig. |H1 
From the hydrodynamical point of view, the less the rj/s, 
the more perfect the matter. In this sense, the nuclear 
matter becomes a more ideal-like liquid around the most 
compressible point in comparison with other evolution 
stages. But note that this minimum rj/s is just a tran- 
sient process. In addition, the extent of approaching an 
ideal-like liquid of the nuclear matter is growing up with 
the increasing of beam energy. In relativistic energy do- 
main, the rj/s of quark-gluon matter becomes very small, 
close to 1/4tt (KSS bound), it is called a kind of perfect 
liquid. 

Temperature dependence of rj/s is an important is- 
sue to understand the transport properties of the nuclear 
matter in different hot and dense environment. To this 
end, we plot a correlation between rj/s and temperature 
in Fig. [Tola) at different energies. Note that the den- 
sity is not fixed in each curve. It is found that there 
is a decrease of rj/s at first when the system is in the 
compression stage. With the system begins to expand, 
rj/s becomes increasing. The higher the beam energy, the 
hotter the nuclear matter, and the smaller the rj/s. From 
this picture it is obvious to find the time when the rj/s 
approaches its transient minimum, is exactly the same as 
the nuclear matter gets the highest temperature. On the 
other hand, in the present beam energy domain below 
200 MeV/nucleon, the transient minimum of rj/s which 
corresponds to the larger compression stage is around 0.4, 
which is about 5 times of KSS bound (i.e. 1/An). 

Furthermore, we can also extract the correlation be- 
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FIG. 9: (Color online) Same as Fig. [5] but for shear viscosity to FIG. 11: (Color online) Same as Fig. [lO]but for density dependence 
entropy density i)/s of the nuclear fireball which is generated in the of rj/s. 
central region. 
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FIG. 10: (Color online) (a): The correlation between rj/s and tem- 
perature at different beam energies with a normal nucleon-nucleon 
cross section parametrization; (b): the rj/s evolves versus temper- 
ature at 130 MeV/nucleon with the different nucleon-nucleon cross 
section. The incident energy and nucleon-nucleon cross section is 
illustrated in the inset. 



tween rj/s and nuclear matter density as shown in Fig. llll 
Here temperature is another hidden variable. Similar to 
Fig. 1101 T)/ s first drops to a minimum value as the den- 
sity is compressed to a maximum point and then rises 
up when the system expands. Larger compressible state 
produces a less rj/s, i.e. the system is close to a more 
ideal-like state. For nucleon-nucleon cross section, the 
effect is mostly reflected before the compressible phase, 
where a larger cross section shows a less rj/s in a given 
density. 

Considering that only final reaction products can be 
detected in experiments, such as the multiplicity and 
flows of the fragments and light particles, it is therefore 
necessary to check the rj/s in the freeze-out stage and 
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FIG. 12: (Color online) The average of rj/s as a function of temper- 
ature at different fixed freeze-out densities (a) and different cross 
sections (b). 



see if it is a useful probe to study the properties of nu- 
clear matter as well as liquid gas phase transition. The 
freeze-out volume has been already studied in many pre- 
vious works [lil [5(1, but in our case, it is more suitable 
to define a freeze-out density instead. The time average 
values of rj/s when the nuclear matter is in some given 
freeze-out density regions of p/po in [0.2, 0.3], [0.25, 0.35] 
and [0.35, 0.45] have been extracted as a function of tem- 
perature. 

Fig. [T2l shows a correlation of the above averaged rj/s 
versus temperature for some given freeze-out densities (a) 
and with different <tnn (b). From Fig. [l"27a) we observe 
that there exhibits a local minimum of rj/s with a value of 
0.7 to 0.8 (about 9-10 times of KSS bound), depending on 
the freeze-out density, in the range of 9 - 15 MeV of tem- 
perature, this phenomenon shall be related to the liquid 
gas phase transition. With the increasing of freeze-out 
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density, we observe the minimal value of 77/s decreases 
and while its corresponding turning or saturated tem- 
perature increases. The former is consistent with the 
results showed in Fig. 12 and the latter can be under- 
stood by the critical temperature/pressure increases with 
the freeze-out density as expected by the pressure-density 
phase diagram [l], 1531 ] . In contrast with the sensitivity of 
77/s to freeze-out density, Fig. [T27b) demonstrates that 
the time averaged 77/s when the system is in a given p/po 
in [0.2,0.4] is not so much dependent on different unni 
which is also consistent with Fig.QTTb). It is interesting 
to note that phase transition temperature 9 — 15 MeV 
which corresponds a local minimum of 77/s is basically 
coincident with previous works (EH, [55| . 

IV. SUMMARY AND OUTLOOK 

Thcrmodynamical and transport properties of a fire- 
ball formed in head-on Au+Au collisions are investi- 
gated in a framework of quantum molecular dynamics 
model. The relationships between different thermody- 
namic quantities are explored. The influences of nucleon- 
nucleon cross section and symmetry energy on the ther- 
modynamical and transport properties are also focused. 
We found that all the properties are very sensitive to 
the nucleon-nucleon cross section and insensitive to the 
symmetry energy. In our calculations, the shear viscosity 
is calculated by a parametrization formula developed by 
Danilewicz and entropy density is obtained by a gener- 
alized hot Thomas Fermi formalism. The present work 
gives a time evolution of shear viscosity over entropy den- 
sity ratio of nuclear fireball, which shows that a transient 



minimal 77/s occurs in the largest compression stage. The 
results at different beam energies show that the larger the 
compression, the more perfect the nuclear fireball be- 
haves. In the present beam energy domain below 200 
MeV/nucleon, this transient 77/s approaches to 5 times 
KSS bound. In addition, temperature and density depen- 
dencies of 77/s are also investigated. It is of very interest- 
ing to observe that a local 77/s minimum, which is about 
9-10 times KSS bound, emerges from the temperature de- 
pendence of 77/s at different constant freeze-out densities 
(025 - 0.4 po), which corresponds to a liquid-gas phase 
transition occurring in the intermediate energy heavy- 
ion collisions. From the temperature dependence of 77/s, 
we learn that the phase transition temperature rises up 
with the freeze-out density. However, the present work is 
still in a phcnomcnological level for investigating 77/s of 
hot nuclear matter which is formed in intermediate en- 
ergy heavy-ion collisions, experimental measurements of 
77/s are still not available so far. Therefore, proposals for 
direct probes of shear viscosity and entropy density in in- 
termediate energy HIC are very crucial and welcome for 
constraining the transport properties of nuclear matter 
around the liquid-gas phase transition. 
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